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Abstract 

We discuss the improvement of flavour non-singlet point and one-link lattice quark 
operators, which describe the quark currents and the first moment of the DIS struc- 
ture functions respectively. Suitable bases of improved operators are given, and the 
corresponding renormalisation factors and improvement coefficients are calculated 
in one-loop lattice perturbation theory, using the Sheikholeslami-Wohlert (clover) 
action. To this order we achieve off-shell improvement by eliminating the effect 
of contact terms. We use massive fermions, and our calculations are done keeping 
all terms up to first order in the lattice spacing, for arbitrary rr? j]?, in a gen- 
eral covariant gauge. We also compare clover fermions with fermions satisfying the 
Ginsparg- Wilson relation, and show how to remove 0(a) effects off-shell in this 
case too, and how this is in many aspects simpler than for clover fermions. Finally, 
tadpole improvement is also considered. 
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1 Introduction 



There has been considerable progress in obtaining reahstic results from nu- 
merical simulations in lattice QCD. A new generation of massively parallel 
computers promises results that can be compared to a wide class of exper- 
imental data. Nevertheless, the finiteness of the lattice spacing always leads 
to systematic errors in the simulations. Therefore there is great interest in 
improving lattice QCD calculations. A systematic improvement scheme, re- 
moving discretisation errors order by order in the lattice spacing a, has been 
suggested by Symanzik [1], and developed by Liischer and Weisz [2] for on-shell 
quantities. 

An 0(a) improved fermionic action which is widely used in lattice Monte Carlo 
simulations is that proposed by Sheikholeslami and Wohlert [3]: 



= E I - ^(^)^(^) - - E ^(^ + «A) ulix) [1 + 7.] ^{x) 

ip{x - afi) U^{x - afi) [1 - 7^] tl){x) 

-2k cs^g ^ ^^(x) a,,Ffr'{x)i;{x] 

where k is known as the hopping parameter and -F^Jf^^"^ denotes the standard 
"clover-leaf form of the lattice field strength.^ If the parameter Cgwi which 
gives the strength of the higher- dimensional operator, is correctly chosen this 
action has no 0(a) errors for on-shell quantities such as hadron masses. For 
perturbative calculations it is simpler to use a slightly different normalisation: 



X 



ij{x - afi) U^{x - afi) [1 - 7^] i){x) 



The parameters of the two forms of the action are related by 



We will use the convention afj,^, = (i/2)(7^7i^ — ju^fi)- 
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where Kc is the critical value of the hopping parameter, at which chiral sym- 
metry is approximately restored. 

Simply improving the action does not remove 0{a) errors from operator matrix 
elements. To do this the operators must also be improved by adding higher 
dimensional irrelevant operators with appropriate improvement coefficients. 
The operators also need to be renormalised. In this paper we will discuss the 
perturbative renormalisation and improvement of bilinear quark operators. 

However, the action (2) with its single tunable improvement parameter Csw 
only improves on-shell quantities. Off-shell quantities still have 0(a) errors, 
which arise from short-distance "contact" terms. We will show how the contact 
terms can be removed at the one-loop level of lattice perturbation theory, and 
off-shell quantities free of 0{a) discretisation errors can be extracted from 
Green's functions. 

There are several reasons why it would be desirable to understand the improve- 
ment of off-shell quantities. In particular the non-perturbative renormalisation 
suggested in [4] involves comparing lattice measurements of off-shell Green's 
functions with continuum perturbation theory results [5] in order to relate 
lattice quantities to conventional renormalisation schemes such as MS. This 
matching will work best at large virtualities, where the running coupling con- 
stant is small, and the effects of non-perturbative phenomena such as chiral 
symmetry breaking have died away. It is obviously desirable to remove the dis- 
cretisation errors in the off-shell lattice Green's functions before making the 
comparison with the continuum. Even within perturbation theory it is easier 
to calculate Green's functions at » than in the region where and 
are comparable. 

Our strategy is to look at the tree-level results for the Green's functions, and 
see what 0(a) effects are present, and what has to be done to remove them. 
We then look at the one-loop perturbative results, and see whether the tree- 
level procedure still works. We find that at one particular value of the clover 
coupling the 0(a) effects arc of the same form as in tree-level, and that then 
we can remove 0(a) effects completely, and find improved Green's functions 
that are free of 0(a) discretisation errors, both on-shell and off-shell. 

Our aim is to find perturbative expressions at one-loop for the MS-scheme 
renormalisation factors and for the improvement coefficients. To do this we 
have to compute each Feynman diagram including all 0(a) terms. These re- 
sults are applicable to both quenched and dynamical calculations of flavour 
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non-singlet matrix elements. 



In this paper we consider the complete set of point operators 



(5) 



with 



= 1, 75, inib, (ynui^- 



(6) 



For one- link operators we discuss the physically interesting case of the leading- 
twist operators occurring in the operator product expansion for the moments 
of the hadronic structure functions [6] . We consider the operators which mea- 
sure the first moment of the unpolarised and polarised structure functions: 



where symmetrisation over ji and v and removal of trace terms is always to 
be understood. 

The perturbative renormalisation of improved point operators has been dis- 
cussed by several groups [7-9]. They use the tree-level values for the operator 
improvement coefficients, q, (defined below) and for the coefficient Csw in the 
Sheikholeslami-Wohlert action. The same settings have been used to calculate 
the renormalisation factors for the one-link [10] and two-link [11] quark op- 
erators in the chiral limit, performing on quark operators the transformation 
discussed in [9]. In this paper we present the Z factors with coefficients q and 
Csw kept arbitrary. This allows us to determine the perturbative contributions 
of the various terms and their relative magnitudes. Moreover, this will enable 
us to implement tadpole improved perturbation theory. 

This paper is organised as follows. In Sect. 2 we give the operator bases for 
the improvement of the point and one-link operators. In Sect. 3 we present a 
method with which to improve the lattice quark propagator off-shell by tak- 
ing care of contact terms, and in Sect. 4 we extend this procedure to improve 
off-shell quark bilinear operators as well. In Sect. 5 we compare with fermions 
satisfying the Ginsparg- Wilson relation, and show how to remove 0{a) effects 
off-shell in this case too. Finally, in Sect. 6 we apply tadpole improvement 
to our perturbative results, and in Sect. 7 we present our conclusions. The 
(sometimes cumbersome) complete results for renormalisation factors and im- 
provement coefficients are collected in the Appendix. 




(7) 
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2 Bases for improved operators 



In this section we write down a general operator basis for the improvement of 
quark operators. The base operators must have the same symmetry properties 
as the unimproved ones, i.e. their transformations under the hypercubic group 
and charge conjugation are determined by the original operator. 

First we consider the five point operators of eq. (5). Subject to the symmetry 
constraints we find the following improved operators: 



fQsy^P = (v>V>)™'' = (1 + amco)ipip - , (8) 

(O^y""" = ('075'0)™'' = (1 + a m CQ)ijj-f5ijj + ]^ac2d^ (^V'T/.TsV') , (9) 

(Y\ imp / — \ imp , , — 1 - ^ 

J = ) = (1 + a m Coj^'Jij.^ - -acitljD^tl: 

+]^a\C2dx(^(T^x'^ , (10) 

(V'TmTsV^)™'' = (1 + amco)^7M75^ 

1 - « 1 /- \ 
--aicii^a^xlbDxi^ + 2«C2'9^ (i^lbi^j , (H) 

/ TT \ imp / - \ imp - 

[Pi^^j = [ilja^^'y^ipj = (1 + amco)ilja^,y-f5ilj 

+^0101-0 {y,j,D^ - -i„D^ + ^aic2 e^^^Ar dr (^V'7aV') , (12) 

where m is the bare fermion mass and D = D — D is the symmetric covariant 
derivative. We have used the lattice definitions 



Un{x)'^{x + a/i) — C/J(x — ajl)'ijj(x — ajl) 
ip{x + afi)Uji^{x) — ip{x — afi)Un{x — ajl) 



(13) 



The C2 terms in the above equations are irrelevant for forward matrix elements, 
which are all that we consider. Therefore we are left with the expressions in 
Table 1. 



We include the terms proportional to cq so that we can get m-independent 
renormalisation constants and at the same time maintain 0{a) improvement. 

Using the scalar operator as an example, there is an equation of motion that 
says that for on-shell measurements ip pip+rjmipip = 0, where t] is a coefficient 
that depends on g. So we can compensate for changes in ci by making changes 
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Improvement basis 


s 


(1 + amco)iptp — ^aciip 


p 


(1 + amco)V'75V' 


V 


{l + amco)tjjjij,ip - ^aciipDij^'tp 


A 


— _ ^ 

(1 + amco)V'7/*75V' - 20'iciilja^x75D xip 


H 


(1 + amco)V'o-//i.75'0 + ^aiciV' (^Tm^;^ " Ti^^mJ 75^' 



Table 1 

The improvement bases for the scalar (S), pseudoscalar (P), vector (V), axial (A) 
and tensor (H) operators. 

in Co, which would allow us to eliminate one of the improvement terms if we 
were only interested in on-shell quantities. This equation of motion means that 
Co is linear in ci if we parameterisc our operators as in eq. (8). If we use other 
parameterisations, for example (1 + abm){ipip — ^ac'iip pip) [12], we would no 
longer find that b was linear in c'l- 

We also consider the conserved vector current: 



+\aic2dx{i)(Ti,x'il^. (14) 

We know that this should need no improvement for forward matrix elements, 
because the only improvement term, C2, is the coefficient of a total derivative, 
and so has no effect on forward matrix elements. Thus provides a useful 
check of our improvement method. 

Next we consider the one-link operators of eq. (7). Here we choose as a basis 
for the improved unpolarised operator 



{O^u)'^^ = (1 + a m Co) l-ipJi^D^ip + ^^^1^1^ V'ct^a \d„, Dx] ip 
-- a C2ip {d^, D^} + - a i C3 ^ dxi^a^xDui^ ■ 



(15) 



This operator basis is the same for the two possible irreducible representations 
of the lattice hypercubic group to which the original operator may belong: 
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T3 (non-diagonal, ji ^ v) and r} ' (diagonal, ji — v). (Our notation for the 
irreducible representations of the lattice hypercubic group follows [13].) In the 
case of the polarised structure function we find that the improvement terms 
allowed by the hypercubic symmetry are different for the representations rf^ 
(non-diagonal, [i ^ v) and r'f^ (diagonal, n = v). When ji ^ v the improved 
operator has the form 



/ \ imp 1 - ^ 1 - ^2 

\Pnu) = (1 + a m Co) -^^7^75-0^-0 - ^ ia ci ilja^^^^^D^il: (16) 
-^aic2 X] '4^(^tiXlb{Dx,D^^i) + ^ac:id^(^'-i5D^^lj\, 

whereas in the traceless diagonal case one has only one improvement term in 
the forward case, so there is no C2, and thus 



^iii,) =(l + amco)-^7;,75£'/,'0--aici 2^^(t^a75 |^a, ^/.j ^ 

A 

^-^^ac■i^^,[i)-i^D^,1^. (17) 

Here repeated // and v indices are not summed over. We will always construct 
a traceless operator when the indices are equal by using the combination 
|'0|(7/i75-D^ — ^v1bDv)4>-i and a similar one in the unpolarised case. 

The coefficients Cq, . . . , C2 can be appropriately determined using the method 
explained in the following sections so that the desired improvement is achieved. 



3 Improving the quark lattice propagator 



3.1 Method 



Even when the fermion action has been improved for on-shell quantities there 
are still 0{a) effects present in off-shell quantities such as the fermion prop- 
agator at a general Euclidean momentum p. In this section we will discuss 
how to find an improved fermion propagator off-shell. First we will look at the 
tree-level Wilson propagator and show how to remove its off-shell discretisa- 
tion errors, then we will generalise this improvement method to the interacting 
case. 

We are used to writing down expressions for 5"^, the inverse quark propagator. 
In the main 0{a) effect is the addition of the momentum-dependent Wil- 
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son mass term. However it is also instructive to look at the quark propagator 
S itself, rather than the inverse propagator. 

Let us start by looking at the propagator at tree-level. From the expression 
for the inverse propagator 



(^*'''*') ^{p,m) = -Y^ 7m sm(ap^) + m + - ^ (1 - cos(ap^)) 



we derive that 



i^+m+ iap2 + 0(a^) (18) 
(i|( + m — ^am'^){l + am){l — ^0(1/ + m)) + O(a^) 



S'-{p, m) = ' "7 ^ + ^ + O(a^) 
ip + m — ^am'^ 2 

= il-am)ST''{p,m,) + ^ + 0{a^). (19) 

The tree-level lattice propagator consists then of two parts, one part is pro- 
portional to a normal continuum propagator with a mass = m — ^am"^, 
and the other part is a momentum independent term. The nature of these two 
parts becomes even more clear when we write them in position space:0 

S'^'^x, y, m) = (1 - am)Sl''%x, y, m^) + ^a5(x - y) + 0{a^). (20) 

We see that (except at short distances, where an additional contact term 
appears) the lattice Wilson- fermion propagator is proportional to S*^, which 
has the form of a continuum propagator with an "improved" mass m^, and 
which has no 0(a) discretisation errors. We will always use -k to mark bare 
quantities which have been 0{a) improved. 

This concentration of the 0{a) effects at short distance is what we should 
expect, in fact the fermion propagator at |a; — ?/| ^ a is an on-shell quantity, 
so it should be automatically improved when the action is improved. It is only 
at short distances of order a that the lattice propagator has a different form 
from the continuum propagator. The necessity of subtracting a 6 function from 
the lattice propagator to obtain an improved propagator has been discussed 
in [4]. 

What should we expect beyond tree level? Let us write the inverse fermion 



^ On the lattice we define 5{x — y) = Sx^y/a'^, where d^^y is the Kronecker delta 
function. 
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propagator as a series in the lattice spacing a: 



S ^(p, m) = cro(p, m) + acri(p, m) + O(a^), (21) 

where the coefficients and (Ti are power series in g"^. On-shell improvement 
tells us that the lattice fcrmion propagator should be proportional to a contin- 
uum fcrmion propagator except at short distances, so wc expect equations of 
the same form as eqs. (19) and (20) to hold, though only at the value of Cgw 
corresponding to on-shell improvement. Thus the propagator should satisfy 



S{p, m) = — — ^ + -A^ + 0{a 

= (1 + ab^m)S^{p, + + 0{a^), (22) 

where the improved bare mass is related to m through 

rrii, = m{l + abmrn) . (23) 

The improvement coefficients 6^, bm and A^ are independent of p and m, and 
should only depend on the coupling constant g'^. By comparing with eq. (19) 
we see that the tree-level values are 

b^ = -l, bm^-1/2, X^^l. (24) 

The propagator S^, is free of 0(a) effects so we call it the improved fermion 
propagator. Later, when we come to define renormalisation constants, we will 
always define them in terms of the improved bare propagator S^, and improved 
bare mass m^. 

Taking the inverse of eq. (22) gives 



5" ^(p,m)^ao(p,rn^)(l-ab^rn)(l-laX^ao(p,rn^fj +0(0^) (25) 
= ^ao{p, m) + abmrri^-^aoip, m)j 

X (1 - ab^m) (1 - ^aX^ao{p, mj^ + O(a^), 

so that dropping terms of order we get 



5" {p,m) = ao{p,m) +a 



|A^ (cTo(p, m))^ - b^mao{p, m) (26) 



+bmm'^-^ao{p,m) 



+ 0{a'). 
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Comparing with eq. (21) we see that the improvement prescription (22) can 
only work if the non-linear relation 

d 

ai{p, m) = -|Ay, ((To(p, m)f - b^mao{p, m) + 6^m^— (To(p, m) (27) 

is satisfied. In subsection 3.2 we shall see that this is indeed the case in one- 
loop perturbation theory. 

The pole mass of the fermion is the p value where CoiP^ vanishes, so for an 
on-shell fermion the factor (1 — |aA^(To(p, m^)) simply reduces to 1. Therefore 
we can see from eq. (25) that the improvement coefficient only has an effect 
when the fermion is off-shell, so we only need to know A^ if we are interested 
in extracting numbers from off-shell lattice measurements. 

From eq. (22) we can find an explicit expression for S^,: 

S*{p,m^)= ^ (s{p,m) -^X^) = V (2^) 

l + amb^\ 2 / (7o(p, m*j 

and 



S^^{p,m^)^{l + amb^) i^S ^{p,m) + S ^{p,m)-X^S ^{p,m)j 

= (7o(p,m*). (29) 

The reason we are interested in S^{pjm^) is that it is a quantity free of 0{a) 
effects which can be constructed from the quark propagator S{p,m), and the 
latter is something we can measure from non-perturbative simulations on the 
lattice. We will find that eq. (27) is satisfied at one particular value of the 
clover coefficient Csw At this Cgw value one can use eq. (28) or equivalently 
eq. (29) to extract (Jq from lattice measurements. The clover action docs not 
have enough tunable parameters to make the off-shell fermion propagator free 
of 0{a) effects, but this does not really matter, because equations such as 
eq. (29) never-the-less allow us to recover the improved off-shell propagator 
from quantities which we can measure. 

Up till now, we have only discussed the improvement of the fermion propaga- 
tor. The propagator and mass still have to be renormalised. The renormalised 
improved quark mass and propagator are given by 



mi?(/i^) 
Sr{p, mR;i/] 



Zm{lJ^) m(l ambm), 



S^{p, m^) 

Zii/J,'^) Z2{iJ?){l + amb^) 



(30) 



1 



5(p,m)-|A^^. (31) 
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3.2 One-loop results for the quark propagator 



We now want to see if the propagator improvement scheme suggested in 
eq. (22) holds in one-loop perturbation theory, and to calculate the improve- 
ment coefficients and renormahsation factors to 0(5'^)- 

Our calculations are carried out in a general covariant gauge, where the gluon 
propagator is 

G(A;) = fc-(l-a)^, (32) 

with k^= sm{ak^/2). The Feynman gauge corresponds to a = 1, the Landau 
gauge to q; = 0. 

We can write the inverse propagator in the form 

S~^{p, m) = m + ly + - (po(p, rn) + api(p, m)^ + 0(a^ g^), 

(33) 

where Cp = {N^ — i)/{2Nc) for gauge group SU{Nc) . Comparing eq. (33) 
with eq. (21) we see that 



(7o(p,m)=iy+m- — -| po(p,m) + 0(g'^), 

lb TT^ 



1 2 9'C 



^i(p, m) = - ^ p,{p, m) + 0{g% (34) 



We also expand the improvement coefficients to first order in g 



2. 



Ib7r2 



b,i, = 



1 ■ 



1 + 



Ibvr^ 
167r2 



0(9'), 
db, + 0{g') 



4. + 0(/) 



(35) 



If we now substitute (34) and (35) into the quadratic equation (27), we find 
that the p functions must obey the following linear condition if the improve- 
ment procedure suggested in eq. (22) works: 
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pi{p, m) + i/po(p, m) + — T^Po(p, 



2 5m 



(36) 



^^A^y + (dx^ -db^)i^ + (dx^ - 2^^ + d^^) 



m 



The explicit expressions for po and pi can be read from the one-loop expression 
for the fermion propagator up to 0(a): 



5* ^(p, m) = i ^+ m + -p^ 



167r2 

,2 



16.64441 - a - 2.24887 c-jw - 1.39727 cl^ + a L(ap, am) 



+„^(i_r(pV 



m 



— m 



'I6 7r2 



11.06803 - 2 a - 9.98679 Csw - 0.01689 c;^^ 
+ (3 + a) L{ap, am) + (3 + a) T {^'^ /rn?^ 

7.13891 - 0.07187 a + 0.48567 Csw - 0.08173 



— ap 



2 9^Cf 



167r2 

— 1(3 — 2q; — 3 Csw) L{ap, am) 



sw 



lamj 



. g'Cp 

16 



6.34664 + 0.14375 a - 1.48503 c^w + 1.28605 ci 



1(3 + 2 a + 3 Csw) L{ap, am) — {a + 3 Csw) T {p^ /m^^ 



.i(3 + 4a-3c,^) — (l-r(pV 



m 



(37) 



am 



2 9^Cf 

16 TT^ 



14.03413 + 1.07187 a + 15.48574 c.^ - 1.52344 cl. 



1 (9 + q; - 6 Csw) L{ap, am) - |(12 + 5q; - Scsw) T [p^/m^^ 



They are 



Po{p,m) = i^ 



16.64441 -a- 2.24887 Cs^^ - 1.39727 c^^ 
+aL{ap, am) + {l — T {t? jmf^ 



+ m 



11.06803 - 2a - 9.98679 c^vy - 0.01689 c 



,2 

sw 



+(3 + a) L(ap, am) + (3 + a) T {fjmf) 



(38) 
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and 



pi{p,m) 



+i m2 



r 7.13891 - 0.07187a + 0.48567 Cs^- - 0.08173 4^ 
-| (3 — 2 a — 3 Cgw) L{ap, am) 
- 6.34664 + 0.14375 a - 1.48503 Csw + 1.28605 

(3 + 2 a + 3 Csw) L{ap, am) - (a + 3 Cgw) T [p^ / 



'm 



-i(3 + 4a-3csw)^ (l- t(pV 



m 



+m 



14.03413 + 1.07187 a + 15.48574 Cs^. - 1.52344 c^^ 



I (9 + a — 6 Csw) L{ap, am) 
-i(12 + 5a-3csw) T(pVm2) 



(39) 



where 



T(x) = ln(l + x)/x , 
L{x,y) =^E - Fo + ln{x^ + y^) 



(40) 



with Fq = 4.369225 ■ ■ ■ and 7^5 = 0.577216 • ■ ■ . Previously [16] we calculated 
the fermion propagator in the limit m^ <^ p^, but eqs. (38) and (39) are valid 
for any ratio /p^ (but a^m^ and a^p^ must both be small). 

Despite the complicated form of eqs. (38) and (39) it can be checked that at 
Csw = 1; and only at this value, eq. (36) is satisfied by po and pi, and hence 
eq. (27) is fulfilled. This allows us to fix the improvement coefficients, which 
in a general covariant gauge have the values 



csw = l + 0{g^), 



Aw, = l + 



16 7r2 



(10.91085 - 1.85625 a) + 0(/), 



16 TT^ 



16.39210 + 0(^^ 



bm- 2 



1 + 



16 TT^ 



22.79406 + 0(^^ 



(41) 



Both b coefficients are gauge invariant.[^ 



^ There was a mistake in bm in [16], which has been corrected here. 
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In addition to the propagator, eq. (37), our calculation also gives us a one- loop 
expression for the critical coupling [14,15] 



1 + 



167r2 



(12.8587 - 3.4333 Csw - 1-4288 c^^) + 0{g^ 



(42) 



In perturbation theory the quark propagator has a pole in the complex mo- 
mentum plane at — —m^oie- look for a value of where S~^{j),m) in 
eq. (21) has a zero eigenvalue. Using eqs. (38) and (39) gives for the pole mass 



"^poie 1 + T7— ^ [1 + «^( Csw - 1)]6 In(ampoie) 



= m 1 + 



X 1- 



am 



16 TT 

16 7r2 



(43) 



(16.95241 + 7.73792 Csw - 1.38038 4^) 



1 + - — 5 f 12.32005 + 18.70718 Csw - 8.23318 c 
167r^ ^ 



sw 



Note that mpoie is gauge invariant, as it should be [17]. At Cqw ~ 1 the pole 
mass is given by 




^poie ( 1 + ^^^^ 6 In(ampoie) j (44) 

23.30995 j 1^1 - ^ 1^1 + 1^ 22. 79406 



The pole mass becomes a function of m^, with the same value of hm as in 
eq. (41). The unwanted am\\i{am^o\e) term vanishes, and so the logarithm has 
the same coefficient as in the continuum. We will see that this is always the 
case, that when Cgw 7^ 1 the coefficients of the logarithm are changed by an 
amount proportional to am, but at Csw = 1 all logarithmic terms have their 
correct values. 

We define our renormalisation constants Z in two different renormalisation 
schemes, MS, and a momentum subtraction scheme, which we will call MOM. 
In both cases we define the Zs in terms of the improved fermion propagator 

The MS renormalisation constants are defined from 

5.(p, = (//2)5ms (p, m^) = (//^)5ms (p, (z^')^*) , (45) 
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where S'ms is the continuum fermion propagator calculated perturbatively in 
the MS scheme at the scale 



(46) 



-rriT 



167r2 



-a + a In 



167r2 



-4 - 2 a + (3 + a) In f + (3 + a)T [p'/m- 



2 1^2 

MS 



In the MS scheme at the scale ji we find for the renormalisation coefficients 
and as defined in eq. (45): 



1 + 



167r2 



2 a ln(a/x) + 16.64441 - 2.24887 c^v 



-1.39727 c;^ - 3.79201 q; 



167r2 



61n(a//) - 12.95241 



-7.73792 CevK + 1.38038 c- 



sw 



(47) 



(48) 



In MOM we define the Zs at the subtraction scale M through 

Zf^^iW) 



(49) 



when = M^. This implies that 



— 1- 



-Tr 



4Ar, 



-Tr 



^S^^{p, m*) 



p^ 



: m 



^MOM(^2) 



(50) 
(51) 



The advantage of the MOM scheme is that all the quantities involved can be 
calculated on the lattice, so it can be used non-perturbatively too. This is 
different from the MS scheme, where we need to compare with a continuum 
quantity which we can only find perturbatively. 



The Zs in the MOM scheme are not simple as in the MS scheme, they still 
have mass and gauge dependences which cancel in the MS case: 
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^•MOM 



7MOM 



(M^) = 1 + 



(M^) = 1 + 



167r2 



16.64441 -a- 2.24887 c^^k - 1.39727 c' 



(52) 



+ 



a L{aM, am*) + a ^ (l - T (mV^^)) 



+ 0(a), 



„ 5.57638 + a + 7.73792 csw- 1.38038 (53) 

16 TT L 

-3 L(aM, am*) - (3 + a) T (M^m^) 

Note however that Z^'^^ becomes gauge independent when the fermion is 

2 

pole" 



on-shell, i.e. at the point = —m? 



The dependence on the lattice spacing a and clover coefficient Csw is the same 
in the MS and MOM schemes, so that the ratio Z^/Z^'^^ is independent 
of Csw This is as it should be, because the ratio ^ms^^mom jg gjjnpiy the 
conversion factor between the two schemes, which can be calculated in the 
continuum, and so should not refer to the lattice in any way. 



4 Renormalisation and improvement of quark bilinear operators 



4.1 Method 



We are interested in calculating the Z factors and improvement coefficients 
for quark bilinear operators. Let us ffi^st set out our notation for a general 
operator. We consider forward matrix elements only, so improvement operators 
proportional to a total derivative will be dropped. 

All the operators in Sect. 2 have the form 

n 

O'^^p = ^po^p + amcQipOip + a ^ CkipQ^'ip, (54) 

k=l 

where tjjOijj is the original unimproved operator, and the ijjQ^ijj are operators 
with the same symmetries as the original, but dimension one higher. Explicit 
expressions for the Q'^ can be found in Table 1. For example, for the scalar 
operator is — p/2. 

We define the flavour non-singlet Green's function in the usual way: 



x,x',y,y' 
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X ( il){x)il){y) ( O + amcoO + aciQ^ + ■ ■ ■ ) ■ip{y')4){x' 
\ V / v.v' 



y^y I i,,^,A 



1 



e 



-ip-(x—x') 



X ^ ^M-\0 + amcoO + aciQ^ + ■■■ + aCnQ'')M-^ 
= T (m-\0 + amcoO + ac^Q^ + ■■■ + aCnQ'')M-^) , (55) 



where M is the fermion matrix, and denotes the Fourier transform from 
position to momentum space. The cq dependence of G^"^^ is simple. We can 
write 



G^""\p,m]Co,Ci, - ■ ■ ,c„) = G'^'"''(p,m;0,Ci,--- , c„) 

+amcoG^""''(p,m;0,0,--- ,0), (56) 

so we only need to give expressions for G for the case Cq = 0. 

Just as we have contact terms in the fermion propagator, we should expect to 
see 0{a) contact terms arising in eq. (55) when x = y or x' = y'. These will 
give rise to a "contact Green's function", C^, which will have to be subtracted 
from the operator Green's function, just as we subtracted a 6 function from 
the fermion propagator. This contact term is thus given by 



- x,x'/ ^ 

(m-^O + OM-^) . (57) 



Since the coefficient of will be 0(a), we only need to calculate it for 
the unimproved operator ipOip, and we only need the leading order in a. The 
Feynman diagrams needed to calculate to 0{g^) are shown in Fig. 1. There 
is no extra calculation involved. All the graphs needed already occur in the 
perturbative expansion of the operator and propagator.^ 

Finding the contact Green's function is simple when we consider point oper- 
ators of the type i/jTiijj where Fj is any 4x4 matrix. Because there are no 
covariant derivatives in the operator, it is unaffected by the averaging over 
gauge fields, and eq. (57) simplifies to give 

C^'{p,m) =riS{p,'m) + S{p,m)ri = {Ti, S{p,m)} . (58) 



A complete listing of the graphs can be found for example in [10], [11] or [18]. 
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Fig. 1. The one-loop lattice Feynman diagrams needed for the contact Green's 
function Co {p, m) , defined in eq. (57) . 

The contact Green's function are more comphcated in the general case. Their 
expressions for the one-link operators are given in the Appendix (Sects. (A. 6) 
to (A.9)). 

So, by subtracting a contact term proportional to C'^(p, m) in a Green's func- 
tion and by choosing appropriately the improvement coefficients, an improved 
Green's function can be obtained. The resulting expression for a renormalised 
improved off-shell Green's function is: 



Zo{ii^; Ci, • • • , c„) 
Z2{iJi?) (1 + amh^) 



(59) 



X 



G (p,m;co,Ci, 



,Cn) - -AoC (p,m) 



The factor l/[Z2(/i^)(l + amb^)] accounts for the wave-function renormalisa- 
tion. The Green's function in eq. (55) depends linearly on the Cfc coefficients, 
while the renormalised Green's function is independent of the c^s. From this 
we can deduce that l/Zo{lJ^] Ci, ■ ■ ■ , c„) must depend linearly on too, so we 
can write 



Zo{p ;ci, • • • ,c„) = 



^0(/^^o,--- ,0) 



(60) 
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where the (s arc coefficients depending on g'^. At the one-loop level, using the 
fact that Z — 1 + 0{g^) and that the (s are Oi^g"^), we can write 

n 

Zo{l^^\ ci, • • • , cn) = ZoiiJ,^; 0, • • • , 0) - ^ CfeCfe + 0{g^) . (61) 

k=l 

Our final formula for the renormalised and improved Green's function is 



1 



Zo(/i^O,-- - ,0) 



Z2(//2) (1 + amb^ 



1 + Efe CfcCfc 



(p,m;co,ci, • • • ,c„) - -AoC (p,m) 



(62) 



The improvement coefficients 6, c, A and C, are independent of the renormali- 
sation scheme, while the renormalisation constants Z are in general scheme 
dependent. Therefore it can be useful to spht the renormahsation and im- 
provement into two stages, and to define "improved bare" Green's functions 

by 



G?(p,mO 



1 



1 



l + amh^ 1 + EfcCfcCfe 

X G^^'^^'ip, m; Co, ci, • • • , c„) - ^AoC^(p, m) 



(63) 



where = m(l + ambm), using the same value for bm as found from the 
fermion propagator (eq. (41)). As in the propagator section, we will use the 
suffix T*r to denote bare quantities which have been improved to 0{a). 

The second step is then to renormalise this improved Green's function multi- 
plicatively. 



Zo{pi';0,--- ,0) 



Gf(p,m.). 



(64) 



It is useful to write corresponding equations for amputated Green's functions 
too. We define the amputated Green's function A in the standard way: 

A'^""''(p,m;co,Ci, • • • , c„) = S'"^(p, m)G'*^'"''(p, m; Cq, Ci, • • • , c„)S'"^(p, m), 

(65) 

where S{p^ m) is the full fermion propagator. The amputation of eq. (65) 
removes all fermion self-energy diagrams from the perturbative expansion of 
A. The one-loop Feynman diagrams for A are shown in Fig. 2. 
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Fig. 2. The one-loop lattice Feynman diagrams needed for the amputated Green's 
function (vertex function), A^™^(p, m; cq, ci, • • • , c„), as defined in eq. (65). 



The improved amputated Green's function, A^, is naturally defined by 

Af (p, = S:\p, m,)G^{p, m,)S:\p, m,). (66) 
Prom eq. (64) we obtain the renormalised amputated Green's function: 

Ag(p,mi?;/x') = ^2(/x')^o(/x';0,--- ,0)Af(p,mJ. (67) 

Substituting eq. (29) and eq. (63) and using eq. (65) we find 



s 1 + amb^ 



1 + Efc CfeCfe 



A'^""''(p,m;co,ci, • • • , c^) 



—-AOS' ^{p,m)C^{p,m)S ^{p,m) 

+ ^A^ {S'"^(p,m), A'^""''(p,m;co,ci, • • • ,€„)} 



(68) 



Similarly to what was done in the case of the propagator, one can now expand 
A'^, and in powers of a and thus obtain a non-linear condition analo- 
gous to eq. (27), from which the improvement coefficients can be derived. In 
the case of the local operators ipTiip, the expression eq. (58) for the contact 
term means that we can write the expression for A^ in the simpler form 
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From eq. (68) we can see that the two improvement coefficients and \o 
associated with the contact terms are only needed for off-shell improvement, 
because the inverse propagator S~^{p,m) vanishes on-shell. 

4-2 Results for point quark operators 

We shall now calculate the matrix elements of all point operators up to 0{g^a) 
including the finite terms. This goes beyond the work of Heathe et al. [7], who 
only considered the g^alna terms. Including all 0{g^a) terms will enable us 
to compute the improvement coefficients to 0{g^). 

The calculations are carried out for arbitrary m? /p^, not just for m? <^ p^ 
as in our previous papers [16,18,19]. We give the results for the amputated 
Green's functions. 

In this section we show how the improvement coefficients and renormalisation 
constants are calculated in one particular case, that of the scalar operator 



The Green's functions and results for the other operators are given in the 
Appendix. We consider forward matrix elements, therefore we drop the total 
derivative terms in the improved bases (although in the scalar case this does 
not make any difference), which now all have the form 



The one-loop expression for the amputated scalar Green's function up to 0{a) 
is: 




(70) 



— J_ — <— — 

= VrV' - -ac\ip{p, rji/j + amcQi/jri/j. 



(71) 



A'^(p,m;0, cf) = 1 - ai^cf 



+ ^— ^ - 11.06803 + 2 « + 9.98679 Csw + 0.01689 ^ 



16 tt' 



sw 



+ cf (-19.17181 + 13.80068 Csv^ - 3.53833 f^,,) 



-(3 + a) L(ap, am) - 3 (3 + a) T (p^/rn^) 
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4q; ■ 



167r2 

16 7r2 



- 1 + T (pVm^) 



6.34664 - 0.14375 a + 1.48503 - 1-28605 ct 



+cf (11.65102 - a - 1.41422 cs,^ + 0.78465 4^) 

3 3 
+ {-^ + a + acf + - Csw) L{ap, am) 



+3 {a + 2>Csw)T(f/ 
+ 



-4a 



m 



rn? + 



15 ^ , c 15 
Y + 10 a - 3 q; q - —Csw 



m 

p2 



1 - T (pVw 



(72) 



+ am 



167r2 



31.06826 - 4.14375 a - 33.97148 Csw + 3.04688 c 



sw 



+cf (-6.08204 + 2.20074 CsvK + 1.44647 4^) 



m 



+ (9 + a - 6csw) L(ap,am) - 2 (3 + a) „ _ 
+2 (12 + 5 a - 3 cf - a cf - 3 Csw) T [p^/m^) 



We only need to give the expression for A*^ when cq = 0, since the full expres- 
sion with non-zero Cq can be recovered by using eq. (56). 

To improve the Green's functions, we need to choose the improvement coef- 
ficients so that all 0(a) terms in the expressions eq. (63) or eq. (68) vanish. 
It is not immediately obvious that this will be possible, because there are 
many more 0{a) terms than there are improvement coefficients, and there- 
fore more equations to be satisfied than there are unknowns. For general Csw 
we can not satisfy all the equations, we can only remove all 0{a) effects if 
Csw = 1 + 0{g'^). In this case we can derive perturbative expressions for the 
improvement coefficients. The results are 



1-cf 



^— ^ (22.79406 - 8.45146cf), 

16 Tl^ 



s g^Cp 

167r2 



16 TT^ 

8.90946. 



(16.39210- 10.88629cf), 



(73) 



All three improvement coefficients are gauge invariant. There is one free pa- 
rameter in this system of equations. The improvement coefficient cf can take 
any value, but once it is chosen, the values of the other improvement coeffi- 
cients are fixed. This freedom comes from an equation of motion, which allows 
us to compensate for a change in one of the improvement coefficients by adjust- 
ing the other two coefficients. For example, there is a particularly interesting 
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s 


^ + 167r^ 


(- 


-12.95241 - 7.73792 csv^' + 1.38038 + Glna/x) 


p 


1 , 9^Cf 
^ + 16n^ 


(- 


-22.59544 + 2.24887 Csw - 2.03602 4^ + 6 In a^) 


V 


16 TT^ 


(- 


-20.61780 + 4.74556 Csw + 0.54317 c|^) 


A 


1 , 9^Cf 


(- 


-15.79628 - 0.24783 Csw + 2.25137 c|,^) 


H 


16 TT^ 


(- 


-17.01808 + 3.91333 Cs^y + 1.972304^ - 2 In a/x) 



Table 2 _ 

The renormalisation constants Z^^ for the point operators (with no improvement 
term added, i.e. = 0). 



value of cf where Xs vanishes, which means that the scalar three-point func- 
tion contains no contact terms, and so even off-shell it is simply renormalised 
by a multiplicative factor. This value of cf is 

= 1 + 7^5.50582 + (74) 

The improvement coefficients Cq and Xg are only defined at Cgw — 1; but 
can also be defined for general Csw values. The result is 



= L (-19.17181 + 13.80068 Csw - 3.53833 c^ J. (75) 

16 TT 

All these results are gauge invariant. 



We calculate the continuum Green's functions (needed for the Z^^^ factors) in 
the MS (minimal subtraction) scheme. In this paper we use a totally anticom- 
muting 75, even when d ^ A. For the scalar Green's function the result in the 
MS scheme is 



167r2 



4 + 2a- (3 + a) In 



-3 (3 + a) T {p'/ml^) 



+ 1 



167r2 



-i + r(pVmy 



(76) 
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cr 


s 


16 TT^ 


(- 


-19.17181 + 13.80068 - 3.53833 c|^) 


V 


16 7r2 


(- 


-9.78635 + 3.41640 CsH^ -f 0.88458 4^^) 


A 


g^Cp 

16 TT^ 


(- 


-19.37225 + 10.31673 - 0.88458 c|^) 


H 


16 TT^ 


(- 


-16.24376 + 6.85531 Csw + 0.58972 4^) 



Table 3 

The improvement coefficients Cf- Tliese give tfie renormalisation constants wlien 
is non-zero, see eq. (61). 





^0 


s 


1 


„S , g^Cp 

^1 ^ 16 TT^ 


(22.79406 - 


8.45146 cf) 


V 


1 


V , g^Cp 

+ 16 TT^ 


(18.14912 - 


6.96476 ci) 


A 


1 


, g'^Cp 


(18.02539 - 


13.36028 c^^) 


H 


1 


M 1 g^Cp 

+ 16 TT^ 


(16.47708 


- 12.86471 cf ) 



Table 4 

Tfie improvement coefficients Cq for general at Csw = 1- 



We can now calculate Z'^^ which is defined by 

Zf-Wzfif,'; cf)As{p, m; 0, cf) = Af{p, m^s) . (77) 

The result is 



1 , 9^ 
MS/, ,2. S\ _ ^ ^ 16 7r^ 



12.95241 - 7.73792 Cg^^ + 1.38038 c|^ + 61na/x) 



1 + ^ (-19.17181 + 13.80068csw - 3.53833c|Jcf 

TTCM • ^''^ 
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Ar 


s 

V 
A 
H 


1 cf + (16.39210 10.88629 cf) 
1 + (13.71395 2.29988 
1 c^* + (11.46014 8.37649 c^^) 
1 - cf + (10.56742 - 5.51435 cf ) 



Table 5 

The improvement coefficients Ar for general c\ at Csw = 1- 









s 
p 

V 
A 
H 


fg^ 8.83678 
-fgi^ 5.21443 
4% 0.22971 

16 TT^ 

4% 1.58146 

16 TT^ 

2^ 1.44071 

16 TT'' 


1 + 5.50582 

1 + fg^ 11.41407 
1 + fg^ 3.08365 
1 + 5.05307 



Tabic 6 

The improvement coefficients Cg and c\ when the contact term Ar is chosen to be 
zero. 



In the MOM scheme we define the Z for an operator O by 



o^(M')Z^°^(M2)Tr f(Ag°'-'^(p))tAo(p)l = Tr [(A© '""(p))^© """(p) 



(79) 



when = M^, where A^"^^ is the operator's Born term. Applying this defi- 
nition to the scalar operator gives 



Tr[A^(p,m;0,0)] 



(80) 



at = M^, so 
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7MOM 



(M2;0) = 1 + 



167r2 



- 5.57638 - 7.73792 cc^ + 1.38038 c' 



-a + 3L {aM, am^) + 3 (3 + a) T (M^m^) 



(81) 



The same procedure can be repeated for all the local operators. In Tables 2-6 
we give the improvement coefficients and MS renormalisation constants for all 
point operators. The MS renormalisation constants are defined by equations 
analogous to eq. (77). The lattice As and MS-scheme As are all given in the Ap- 
pendix. In several of the tables there is no entry for the pseudoscalar operator 
- this is because it has no Ci improvement term, so the associated quantities 
are not defined. The Zs in the MOM scheme are given in the Appendix. 

The Cq values in table 4 agree with the values given in [20] for the case cf = 0. 



4-3 Results for one-link quark operators 



Wc consider now the operators in eq. (7). We study the improvement of these 
operators along the same lines used for the point operators. The expressions 
for the contact Green's functions (57) will be more complicated, and are given 
in the Appendix. 

From eqs. (15), (16) and (17), we can see that in forward matrix elements 
a basis for the improvement is given in the unpolarised case {0^1,^7^^ and 



(Cm^^)™^ = ^ (1 a m Co) ip-f^^D^ip + ^ a ici ^ ipa^^x [d^, D^] ip 

A 



(82) 



in the polarised case with ji^ v (0L, rf^) by 



/ c \ imp 1 - ^ 1 - ^2 

[Pnu) = -{1 + amco) t/'7^75L'i.V' - ^ ia ci i/ja^^D^ilj 



and in the polarised case when ji — v (OL; '^a'^) by 
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(t- \ imp 1 — <-» 1 — r^*^1 

=-(l + amco)^7M75^:.^- gaici^/^a/,A75|^A,^v|^. (84) 





Z^(/x2;0,0) 


05 ^(6) 
05 ,(3) 


1 + (-1.27959 + 3.87297 Csw + 0.67826 - f Ina/x) 
1 + (-2.56184 + 3.96980 Csvv- + 1.03973 4^ - f Ina/x) 
1 + ff^ ( 0.34512 + 1.35931 Csw + 1-89255 \^ Ina/i) 
1 + ( 0.16738 + 1.24953 Csw + 1-99804 c|^^. Ina/x) 



Table 7 

The MS-scheme renormalisation constants for the one-hnk operators. 





Ci 








fig^ (-4.27417 + 1.08793 c^^^) 








figl (-4.27417 + 1.08793 Csvk) 


05 


^4 


16 TT^ 


(-5.61603 + 4. 10778 cs^y - 0.26315 c|^) 


05 


^(3) 
^4 


16 TT^ V 


-15.31376 + 8.54773 03^^- 0.26036 cl^) 



Table 8 

The improvement coefficients ^i for the one-link operators. 





C2 


^(3) 
05 ^(6) 
05 ,(3) 


( 9.40584 + 4.60327 Csw + 0.46669 4^) 
( 6.67330 + 4.53710 Csw + 0.44621 4^) 
fgl (-8.29791 +4.21724 Csw - 0.49384 4^) 



Table 9 

The improvement coefficients (2 for the one-link operators. 
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Co 


'^^IV 1 1 1 

05 ^(3) 


1 C2 + (16.34500 12.24534C2) 
1 C2 + (17.20377 8.69045 C2) 
1 C2 + (16.28373 10.73103 C2) 
1 - ci + (16.95724 - 11.33434 ci) 



Table 10 

The improvement coefHcients cq for general C2 at Cs^^. = 1 for the one-link operators. 





Cl 


^(6) 

^ixu ) ' 1 
05 ,(6) 

05 ^(3) 


C2 + ^ (1.18321 + 0.84682 C2) 



Table 11 

The improvement coefficients ci for general C2 at Csw = 1 for the one-link operators. 







,{6) 
^/J-vi '3 

^(3) 
05 ,{6) 
05 ,(3) 


1 C2 + (8.13135 3.04784C2) 
1 - C2 -h (9.28735 - 0.52613 C2) 
1 C2 + (8.48845 1.84428 C2) 
1 Cl + (7.96628 3.12321 Cl) 



Table 12 

The improvement coefficients \o for general C2 at Csw = 1 for the one-link opera- 
tors. 

For the one-link operators, we calculate the Green's functions in the limit 
-C P^, keeping terms up to first order in m. Using the results of our 0(a) 
calculations which we have collected in the Appendix, we can derive the values 
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Co 


Cl 


C2 






-16?^ 0-98385 


C2 + 0(5') 


1 + 


g^Cp 

16 TT^ 


5.08351 


o 




0.24789 

16 ■K-' 


C2 + 0(5') 


1 + 


16 TT^ 


8.76121 


^ jjiv-i 


^(6) 


4% 1.09147 

16 TT^ 


1 + fgf 8.67421 


1 + 


16 TT^ 


6.64417 


Qh 


^(3) 


+4% 0.77983 

16 TT^ 


1 + fgf 4.84307 









Table 13 

The values of cq, ci, C2 for which the one- link operators have no contact terms 
(Ao = 0). 

of the renormalisation constants and improvement coefficients which are given 
in the Tables 7-13. For each operator there is a particular value of the qs where 
the coefficient Xq vanishes, and therefore there are no contact terms and even 
off-shell the operator is simply renormalised by a multiplicative factor. These 
values are given in Table 13. 

Note that in the unpolarised case we can only determine Ci to 0{g^) from our 
one-loop calculation. This is because Ci is the coefficient of an operator which 
vanishes at tree- level (because it involves [Di,, D)]). However, we do still know 
the improved Green's function to 0{g^). 



5 Off-shell improvement for Ginsparg- Wilson fermions 

Like clover fermions, Ginsparg- Wilson fermions are free of 0(a) effects on- 
shell. So it is instructive to sec what happens when our off-shell improvement 
conditions (28) and (63) are applied to Ginsparg- Wilson fermions [21]. 

The defining Ginsparg- Wilson relation is 

Daw 75 + 75 Daw = aDow 75 Dqw , (85) 

where Dqw is a Ginsparg- Wilson fermion matrix. Prom this matrix we can 
(at least in principle) define a related matrix [22] 

Kgw ={i- ^Dgw) ' Daw (86) 

The eigenvalues of Dq^ fie on a circle of radius 1/a and centre 1/a, while the 
eigenvalues of Kqw lie on the imaginary axis. Prom eq. (85) and eq. (86) we 
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find that 



Kgw 75 + 75 Kgw = 0. (87) 

The propagator which we would really like to know is the fermion propagator 
corresponding to Kgw- It should have the correct chiral properties, and be 
free from 0{a) discretisation errors. However, we cannot work directly from 
Kgw: because it is non-local. Therefore, we need to write down a formula for 
the propagator we would get from Kgw, but expressed in terms of Dgw- This 
propagator will satisfy chirality even at zero distance, so we expect it to be 
improved off-shell too. 

Let us now add mass to the problem in the same way as it is added in the 
clover case, simply by adding a mass term m.'ipip to the action, giving 

i> [Dgw + m]^ (88) 

as the fermionic part of the action. Another way to add mass effects would be 
to use the alternative action 

iP [(1 - am^/2)DGw + m^] t/-, (89) 

which has the advantage that there is no mass improvement needed. That is 
the method wc used in [21]. Here we have added the simple mass term, eq.(88), 
because we want to preserve the analogy with the clover action. 

If we reexpress the unimproved massive propagator {Dgw + m)~^ in terms of 
Kgw, we find 



. Dgw + m/ a 

= (TT>f( 



Kgw + m/{l + ^m) 
Fourier transforming, we get 



+ 



2 -|- am 



S{x - y). 



'^GW 



(p, m) 



S^{p, m*) + 



2 -|- am 



where 



m 



m^ 



1 + |m 



(90) 



(91) 



(92) 



and S-^. is the Fourier transform of {{Kqw + ^*) ^)a- Eqs. (91) and (92) are 
the analogues of eqs. (22) and (23) respectively, remembering that terms of 
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O(a^) are dropped in Sect. 3.1. Solving eq. (91) for 5"* we find 



S^{p,m^) = (1 + |m)^ 



{p, m) 



2 + am 



(93) 



which has the same form as eq. (28). Note that the only matrix we need to 
invert to calculate this improved propagator is the matrix {Dqw + "^)) which 
is well-defined, and local (in the sense that its elements decrease exponentially 
with separation). Comparing these formulae with those in Sect. 3.1 we see that 

= 1, h^ = -1 and bm = ■ (94) 

These results are independent of g^. These all-order results coincide with the 
tree-level limit of the clover fermion result eq. (41). The values depend on 
the fact that in this paper we have added mass term to the Ginsparg- Wilson 
action in the same way as to the clover action. 

Next we want to improve the Green's function corresponding to a flavour non- 
singlet operator O = ipOip, where O includes Dirac structure and covariant 
derivatives. We want our improved Green's function G^{p, rrii,) to be given by 

Gf(p,m,) = ^ ( \ O ^ ) , (95) 
\ Kgw + Kgw + m^/ a 

where denotes the FoTirier transform. However, we need to re-express it in 
a form that involves only Dgwi not Kgw- This can be shown to be equivalent 
to the expression 



Gf(p,mO = (l + ^m)^ 



(96) 



where 



Gg7(p,m)^^(;^^0--^^) , (97) 

\ Dgw + ^ J^GW -r vTi I A 



= O + amcoO - -c^{DgwO + O Dgw) + -^DqwO Dgw. (98) 

C'^ip^m)^ :F lo ^ + n ^) ' (99) 

\ Dgw + "fn Dqw + /a 



with 



Co = l-Ci, (100) 
Ao = l-ci. (101) 

Eq. (96) has the same form as eq. (63) (up to terms of 0(a^)). A more general 
form of eq. (96) can be found in [21]. 
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Comparing eqs. (100) and (101) with Tables 4 and 5, we see again that the all- 
orders Ginsparg- Wilson result is just the tree-level result for clover fermions. A 
particular point to note is that the Ginsparg- Wilson improvement coefficients 
are the same for all operators, while the clover action improvement coefficients 
are operator-dependent. A further simplification in the Ginsparg- Wilson case 
is that there are no ( coefficients needed, they are all zero. This means that the 
renormalisation constant Zq is independent of the improvement coefficients 
Ci, (see eq.(60)). 



6 Tadpole improvement 

6.1 Analytic results 

It is well known that many results from (naive) lattice perturbation theory 
are in poor agreement with their counterparts determined from Monte Carlo 
calculations. One main reason for this is the appearance of gluon tadpoles, 
which are typical lattice artifacts. They make the bare coupling g into a poor 
expansion parameter. Therefore, it was proposed [23,24] that the perturbative 
series should be rearranged in order to get rid of the numerically large tad- 
pole contributions. This rearrangement will be done by using the variable uq, 
derived from the measured value of the plaquette 

Its value depends on the coupling g'^ — 6/l3 where it has been measured. 

There are two main steps involved in tadpole improvement. 

Firstly, we know [24] that in the mean field approximation the Z for an oper- 
ator with no derivatives is 

Zo ^ Mo""^ , (103) 

so it is reasonable to hope that a perturbative series for (ZoUq°~^) will con- 
verge more rapidly than a series for Z^ itself. Secondly, instead of writing our 
series in terms of the bare parameters, we reexpress it in terms of the tadpole 
improved, TI, parameters 



all 


2 -4 

= 9 , 
















= Ci Uq, 



(104) 
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where n is the difference between the number of covariant derivatives in the 
higher dimensional operator muhiplying Cj and the number of covariant deriva- 
tives in the operator to be improved (n is always 1 for our choice of improve- 
ment terms). The new couphng gf-n is called the "boosted" coupling constant. 
Other choices of boosted couphng are possible, for example one could also use 
the renormalised coupling constant at some scale ~ 1/a^. To carry out this 
rewriting of the series, we simply replace every by [s'ti'^oIs'ti)], where uo{gTi) 
is the perturbative expansion for uq: 

«o(5t:) = 1 - ^^F^' + 0{g'J . (105) 

Formally, this cannot change the all orders result, but it should improve the 

rate with which the series converges. The same procedure is followed with the 
improvement coefficients Csw and q, for example Csw is to be replaced by 

In this paper we will look at the tadpole improvement for operators with no 
anomalous dimension. The interplay between tadpole improvement and the 
renormalisation group, needed when considering operators with an anomalous 
dimension, will be considered in a future paper. The result of this procedure 
is rather simple for the one-loop Z factors. If the original Z is given by 

Zo^l + Bo{ cs^, Ci) + 0{g^) , (106) 

lo vr"^ 

then the tadpole improved Z is given by 



7TI _ l-n£) 



1 + 



9ti 

167r^ 



■Cf {Bo{ cll, cf) + (1 - n^)7r^) + 0{gt,) 







Y _|_ ri DTI I /o^^4 



IGtt' 



,C^55' + 0(y^J 



(107) 



For the V and A operators [ud — 0) in the MS scheme we get the following 
tadpole improved Bq terms: 



^Ti,Ms ^ _iQ_748i9 + 4.74556 cf^ + 0.54317 ( cf^f 

+ (9.78635 - 3.41640 cf^ - 0.88458 ( cf^f) , (108) 

= -5.92668 - 0.24783 c^J^ + 2.25137 ( cf^f 
+ cf (19.37225 - 10.31673 c^i^ + 0.88458 ( cf^f) . (109) 



Tadpole improvement is not just applicable to renormalisation factors - it can 
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also be used to give improved values for the improvement coefficients. The 
improvement coefficients for the fermion propagator, eq. (41), become 



— Mq ^ 



9l 
167r2 



Cf (1.04124 - 1.85625 a) + 0{g. 



9ti 



1 + 6.52250 + 0(^^J 



IGtt 



IGtt^ 



Cf 12.92445 + 0(^^J 



(110) 



An unfortunate ambiguity is that there is of course considerable freedom in 
the choice of boosted g. At one-loop none of the numerical coefficients are 
affected by this choice, so if one prefers another boosted g, all the formulae 
in this section can still be used, the only change is that every g^^ has to be 
replaced by the alternative boosted g. 



6.2 Comparison with non-perturbative results 



To test the validity of tadpole improved perturbation theory, we will now 
compare our results with known non-perturbative results in the quenched 

theory. The local vector current is best suited for this purpose, because the 
renormalisation constant Zy and improvement coefficient Cq are known non- 
perturbatively for a wide range of values of . 

The comparison is done at = 1. At this value of the coupling one finds 
non-perturbatively [25] Csw = 1-769. We will use this number in both the 
perturbative formulae and the numerical calculations. For uq we obtain the 
value 0.8778. We then get 



Z;^' = 0.8242 + 0.04860}", (111) 
4'^' = 1.2733 - 1.0990c[' . (112) 

In Fig. 3 we show and Cq'" as a function of cY ■ We compare the results 
with the numbers of three independent non-perturbative calculations. The 
first calculation [27] uses the nucleon matrix element of the local vector cur- 
rent to determine Zy and Cq . The second one is based on the Schrodinger 
functional [26] , and the last calculation [28] uses chiral Ward identities to im- 
prove the current and renormalisation following [4,5]. In the latter case we 
calculated only at the value of cY where Ay = 0, as in Table 6. It should 
be noted that the results still have errors of O(a^), which can be as large as 
10% [15], so that we cannot expect the results to agree completely. For Zy we 
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Fig. 3. The renormalisation constant Zy and the improvement coefficient Cq as 
a function of c( . The sohd lines are the results of tadpole improved perturbation 
theory. The dashed lines refer to the non-perturbative results of [27], and the symbols 
mark the non-perturbative results of [26] (solid circle) and [28] (solid diamond). 



find good agreement between the tadpole improved perturbative numbers and 
all the non-perturbative results. For our numbers agree with [27] and [28]. 

The Schrodingcr functional result, on the other hand, lies 10% above the 
other numbers. (It is important to remember that different definitions of Z 
may give results differing by O(a^), so both results could be consistent.) 

In Fig. 4 we show the renormahsation constant Zy as a function of g^. At 
smaller values of the coupling (higher values of 0) the agreement between 
tadpole improved perturbation theory and non-perturbative results becomes 
even closer, as one might expect. In those cases where we could check this, we 
found the discrepancy to reduce to ~ 4% at /5 = 6.4. Thus we may say that the 
non-perturbative results agree with those of tadpole improved perturbation 
theory within the expected O(a^) and 0{g^') uncertainties. 
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g 

Fig. 4. The renormalisation constant Zy as a function of . The solid hne is 
the result of tadpole improved perturbation theory, while the dotted line shows 
the one-loop perturbative result with no improvement. The dashed line is the 
non-perturbative result of [26]. 

7 Conclusion 

In this paper we have presented extensive one-loop perturbative calculations 
of lattice Green's functions, in which we have kept all 0{a) terms. This allows 
us to investigate operator improvement, firstly to see what sort of improve- 
ment terms are needed, and secondly to calculate values of the improvement 
coefficients. 

We find that we can produce off-shell 0(a) improved Green's functions, to all 
orders in the Ginsparg- Wilson case, and at least to Oi^g^^ in the clover case. In 
our one-loop calculations we find that we only need gauge-invariant improve- 
ment terms. No extra improvement terms associated with BRST symmetry 
are required at this level. 

Off-shell improvement doesn't mean that there are no contact terms. As long 

as we know the form of the contact terms, we can remove them by using the im- 
provement coefficients A. Contact terms, responsible for the off-shellness of the 
propagators and Green's functions, can be removed using a well-determined 
procedure. There are always particular values of the improvement coefficients 
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for which the contact terms vanish, so that one still has a multiplicative renor- 
malisation. 

In the Ginsparg- Wilson case improvement is particularly simple, because the 
improvement coefficients are universal, they do not depend on the operator 
considered, the coupling constant, or even on which theory we are simulating 
(we assume that the bosonic sector has no 0(a) discretisation errors). This 
is not so in the clover case, the coefficients depend on the coupling, and are 
different for each operator. 

We have the tadpole improved one-loop values for Z factors calculated at 
arbitrary c^n-, and for improvement coefficients calculated at c^w = ^ + 0{g'^), 
which is the only place where 0{a) improvement is possible. Numerical test 
cases show that tadpole improvement works well down to P ^ 6.0 for operators 
with no anomalous dimension. We are investigating tadpole improvement in 
the case of operators with an anomalous dimension. 
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A Appendix 

In this Appendix we give the perturbative expressions for the amputated three- 
point functions (vertex functions) for all the operators we have considered 

(apart from the scalar density, which is given in Sect. 4.2 of the main text), 
calculated to 0(a), and the values of their improvement coefficients. In order 
to make transparent the transformation of these numbers into the popular MS 
scheme the corresponding continuum quantities are also given. 

In order to shorten the expressions for the Green's functions we will use the 
functions T and L defined in eq. (40). 
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A.l Pseudoscalar Vertex 



The pseudoscalar operator is simpler than the other operators because there is 
no Ci improvement term possible (and also none needed) for the forward three- 
point function, because p and 75 anti-commute. The one-loop expression for 
this vertex up to 0(a) is: 



A^(p, m; Co = 0) = 75 75 



16 7r2 



- 1.42500 + 2 q; + 3.43328 c 



sw 



— {2>-\- a) L{ap, am) — {3 + a)T (^p^/m^^ 



+am 75 



167r2 



3.82788 - 0.14375 a - 2.49670 ci 



+(3 + a)L{ap, am) + 2 (3 + a) T (^/m^) 
In the MS scheme we have 



(A.l) 



^Ms) =75 + 75 



167r2 



4 + 2 a - (3 + a) In j 
-(3 + a)r(pVmy 



(A.2) 



The MOM scheme renormalisation factor is defined by 



7M0M 



(M^)Z^o^(M^ 



IV[75A^(p,m;0)] 



at p^ = M^. The result is 



(A.3) 



7MOM 



(M2) = 1 + 



g'c, 



- 15.21941 -a + 2.24887 Csw - 2.03602 c 



167r2 

-|-3 L {aM, am 



SW 



(3 + a) T (mV 



+ 0(a). 



(A.4) 



As well as the lack of a derivative improvement term, another special fea- 
ture of the pseudoscalar operator is that the improvement terms amA^ and 
a{S~^, 75}, which appear in eq. (69), have the same functional form to this or- 
der, so there is no natural way of determining Cq and Ap separately. We choose 
to improve the operator by setting Ap = 0, and making all the improvement 
through the mass- dependent term. Then the improvement coefficients are (at 

C-sw — ^) 
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9^ 
167r2 

Ap = 0. (A.5) 



Co =-f^ 5.21443, 



A. 2 Local Vector 



The one-loop expression for the local vector vertex up to 0{a) is: 



Ajf (p, m; 0, c[) = 7^ - a ip^ c\ 



V 



+ 7m 



167r2 



3.97338 + a - 2.49670 c^^ + 0.85410 c^^ - (9.78635 



-3.41640 Csw - 0.88458 4^) - a L{ap, am)-a^(l-T [p^/m^)) 



+ 



+ 1 



16 TT^ 



-2a + 4a^(l- T(pVm2)) 



p 



+ aip^ 



P^ 16 TT^ 



2 (3 + a) (l - T (pVm^ 



16 7r2 



- 8.66505 + 2.85625 a + 9.52789 Csw - 0.39053 c 



sw 



+c[ (18.59361 -a- 2.24887 - 0.16098 c; 



swJ 



{3 — a — cY a — 3 Csw)L{ap, am) + 2(3 + a) 



m 



+ p^ 



-(18 + 8a-6cY -SacY -3 Csw) ^ (l - T (pV^^)) 



16 7r2 



7.64168 + 0.85625 « + 7.74287 - 1.38589 



swJ 



+c[ (13.96523 ~a- 5.54358 c^^ - 0.36162 c; 
— ^ (3 — 6 cY — 2 a — 3 Csw) L{ap, am) 

-(3-3csw-^cY -acY)T (f/m^) 

2 

+1 (3 + 6a-3csv^-2acl') ^{l-T{f/m^)) 



(A.6) 



3 + 4Q;-2c]:'Q; + 3csv^-2(a + 3csw-)T(pVm2) 
p^ Ibn^ I ^ ^ 



m"' + p"' 



— I 6 + 12 q; — 4 c^ q; — 6 Csv 



p^ 



(l-T(pVm^)) 



In the MS scheme we have 
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167r2 



a — ahi 

2 



P 



+ 



+ 1 



16 TT^ 



2Q; + 4a^ 



l-T(pVm|3)) 



2 (3 + a) (l-r(pVmy) 



16 TT^ 

We consider two ways to define Z^'^^ for tlie vector, 



(A.7) 



1 



1 



AN, 



PjiPv 

y^^Tr 



■Tr b,K] 



1 



7M0M 7M0M ' 

vtrans 



7MOM yMOM ' 

^2 ^Mong 



SO tfiat 



(A.8) 
(A.9) 



= 1 + 



167r2 



(A.IO) 



«M^;0) = 1 + 



20.61780 + 4.74556 c^w + 0.54317 cl^ + 0(a), 
'f 



167r2 



-20.61780 + 2a + 4.74556c 



+0.54317 _ 4a ^ (1 - T (MVm^ 



(A.ll) 
+ 0(a). 



A. 3 Conserved Vector 



The one-loop expression for the conserved vector vertex up to 0{a) is: 



+ 7^ 



16.64441 + a + 2.24887 c™ + 1.39727 c^. 



im 



+ 



+ 1 



167r2 

a L{ap, am) — a — (l — T [p^ / 

-2a + 4a^(l-T(pV 
2 (3 + a) (1- r(pV^')) 



sw 



p"^ 16 TT^ 



m 



P^ 16 TT^ 
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167r2 



11.27782 + 1.85625 a + 3.97134 c.w - 0.16345 c 



m 



— (3 — 2a — 3csw)L(ap,am) + 2(3 + a) „ „ 

jj^Z _j_ pZ 



-(12 + 5 a - 3 cs^) ^ (l - r 



167r2 



6.34664 - 0.14375 a + 1.48503 Csv^- 1.28605 



+1(3 + 2 q; + 3 Csw) L(ap, am) + (3 Cgw + 0()T (p^/rn^^ 
+ 1 (3 + 4 a - 3 cs^) %{l-T (pVm^)) 



(A.12) 



+ a 



P^ 16 TT^ 

i2 



3 + 2q; + 3 Csw — 2(q; + 3csw) T [j? jvr?^ 



+4 a 



6 + 8Q;-6cs,^j — (l-r(pV 



In the MS scheme we have 



167r2 



q; — q; In 

2 



-a^(i-r(pV-y) 



+ 



+ 1 



P^ 16 TT^ 
P^ 16 TT^ 



2a + 4a!5^(l-T(pVmy) 



2 (3 + a) (l-r(pVmy) 



(A.13) 



We can define longitudinal and transverse Z^'^^ according to eqs. (A.^ 
and (A. 9), giving us 



«M^) = l + 0(a) 



g'Cp 

167r2 



2a-4a^(l-T(MVm; 



(A. 14) 
+ 0(a).(A.15) 



The conserved vector current satisfies eq. (63) without any improvement terms, 
i.e. 



A. 



0, 
0. 



(A.16) 
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This is as it should be, the conserved current is aheady improved for forward 
matrix elements, and no further improvement terms are needed. 



A. 4 Axial Vector 



The one-loop expression for the axial vector vertex up to 0{a) is: 



A^(P, rn; 0, cf) = 7^75 - a^(^7^75 + 7/.75^)cf 



16 7r2 



- 0.84813 + a + 2.49670 Csw - 0.85410 c 



sw 



-4 (19.37225 - 10.31673 c^v^ + 0.88458 c^.,^,) - a L{ap, am) 



-2(1 + a)T {p'/m') - (2 - a) ^ (l - T {p'/m')) 



p. 



16 TT^ 



-2 a -4(1 -a) T(f/m^) 



+4(2-a) ^(l-T(pVm^)) 



+ 2(^W + 7m75:^)^Y6^ 

+ a^(/7M75 + 7M75/)Y^ 



2 (1 - a) (1 - T (pV^^)) 



1.34275 + 0.85625 a - 1.71809 



+0.13018 cl^ + cf^ (13.96522 - a + 0.54301 c^^^ + 0.05366 c; 
+Q; (1 + ) L[ap, am) 

+2(0; + Csw) T (f/m^) +2(1-0;) 



sw) 



m 



+ p2 



;i - 4a - 2c^^ + ac^^ + 2 c^v^) ^ (l - ^ (pV^')) 



+ am-f^-f5 



16 7r2 



7.47851 - 1.14375a - 8.742870^^. + 1.38589c 



2 



+cf (10.92421 -a- 4.44321 c^v^ + 0.36162 c; 
+1 (^6c^ + 3 + 2q; — 3 Cgvi^j L{ap, am) 

+(3 + 6 a + (1 - a) c^^ - Csw) T i//m^) - 4 



sw) 



m 



m^ + 



+i(l3_6a-2(2-a)cf- c,^.) ^ (l - T (pV 
-3 + 4Q;-2Q;cf+ Cg^^. + 4 



Pm ^ ^^C'f 
+ am ^-75—^ 

lb TT^ 



m 

2 



+p^ 
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+2 (6 - 5 a - 2 (1 - a) - Csw) T (p'/m') 



-2(13-6a-2(2-a)ci^- c 



m 
sw) 5" 

pi 



(l-T(pVm2)) 



(A.17) 



In the MS scheme we have 



167r2 



a — q; In 

2 



+ 



-2 (1 + a) T (//my - (2 - a) ^ (l - T (pV^ls)) 

16 TT^ 



-2a-4(l-a)r(pVmy 



+4 (2 - a) ^ (l - T (pV^Is)) 

+ ^(^7/.75 + 7m75^)^ 



2(1 -a) (l- T(pVm|j3; 



(A.l^ 



Again we can define both transverse and longitudinal Z^'^^, 



si^'V j 4A^/'L^^V^-J-Z2^0M- 



r P' 47V, 



75^^ 



^^■MOM^-MOM • 

-'4.trans 



^'MOM^'MOM" 

2 -^long 



(A.19) 
(A.20) 



This gives the MOM scheme renormalisation factors 



Z^o^(M^O) = 1 + 



^long 



15.79628 - 0.24783 Csw + 2.25137 



16 7r2 

+2(l + a)r(MV7 
+2(l-a)^(l-T(MVm^)) 



+ 0(a), (A.21) 



- 15.79628 - 0.24783 Csw + 2.25137 c 



167r2 

+2 a + 2 (3 - a) T (mV^^) 
-2(3-a)^(l-T(MVm^)) 



+ 0(a). (A.22) 
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A. 5 Tensor vertex 



The one-loop expression for the tensor vertex up to 0{a) is: 



16 7r2 



4.16568 - 1.66446 Cc 



0.57503 c: 



- cf (16.24376 - 6.85531 Csw - 0.58972 c^ J + (1 - a)L{ap, am 
-(1 -a)T (f/m^) + 2(1 - a) ^ (l - T (pV^"')) 

+ i 4( W - 7^P.)75 [2(1 -a)T (f/rn') 



4(i-a)^(i-r(pV-^)) 



4 (1 - T (pV^^)) 



+ a^(/cr^^75 + f^M,.75#)77— I 

Z 16 TT^ 



-0.42868 c^ 



- 3.66115 + 1.85625Q; + 0.96286 Csw 

iO TT" L 

cf (16.27942 - a + 0.26479 Csw - 0.26155 c^,^) 



+1 (2 q; cf — 3 + 2 Q! + Csw) L{ap, am) + {a — Cgw) T [p^ /m^) 

2 

+1 (2 (4 + a) cf - 19 - 4 q; + 3 Csw) "^(l-T (p'/m^)) 

2 

m 

+4- 



+ am(7„^75^^ — ^ - 7.42480+ 1.85625a + 5.16191 c^^^ - 0.09170 c^^ 

16 TT"* L 

+ cf (17.60663 -2a- 7.02465 Csw - 0.24108 c^^) 
— (2 — q; — 4 cf — Csw)L{ap, am) 

+2 ((1 + a) cf - a + cs^)) T (p7m') 

2 

-(5 - 6a - 2(1 - a) cf + c,^) ^ (l - ^ (pV^')) 



+ iam^(7,p^-7^p,)75^ 



3 + 



+2 ((1 - a) cf - (3 - 2q; + Csw)) T (f/m") - 2(1 - a) 



m 



m^ + p^ 



+2 (-2(1 - a) cf + 5 - 60; + Csw) ^ (1 - ^ (pV^')) 



(A.23) 
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In the MS scheme we have 



^is(P> ^Ms) = f^Mi'Ts 



167r2 



In 



+2^(l-T(pV-y) 



+i^(7i^Pm - TmPj^JO's^^^ (1 - a 



9 



16 7r2 
(l - T (//my) 



2T(//my 



16 TT^ p2 



1 - T (pv^y 



(A.24) 



We define Z 



MOM 



by 



12 V 47Ve 



75 (7^^ A^^ 



yMOM yMOM 
^2 ^i/ 



(A.25) 



at the scale — M^, which gives 



^MOM(^2. 0) = 1 + 



16 7r2 



20.81009 + 3.91333 Csw + 1-97230 c;^ 

2 



+ a - L (aM, am*) - a ^ (l - T (MVm^)) 



(A.26) 



v4.5 First unpolarised moment, off-diagonal (ix^ u), symmetrised 



For the one-hnk operators, we calculate the Green's functions in the limit 
<^ P^, keeping terms up to first order in m. This is sufficient to calculate 
the improvement coefficients. 

Our improved operator for the first unpolarised moment in the t^"^ represen- 
tation of the hypercubic group is 



(6) 1 - 1 / <-> \ 

0]l =(l + amco)-0-(7^L'. + 7.^M)V' 



2 " 2 
-1 

'2 



1 -1 ^ ^ ^ 
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8 

The result for the amputated operator Green's function is: 



(6) 1 

A's (p, m; 0, ci, C2) = i- (7/,^^. + 7^P/,) + C2ap^p^ 



[--a) InaV - 18.80927 + 3.79201q; 



-1.62411 csw + 0.71900 + ci(-4.27417 + 1.08793 Csw) 
+C2(-9.40584 + 4.60327 Csw + 0.46669 4^) 



+ i 



+ aPuPv 



16 TT^ 

'f 



-m 4 
16 vr^ 



167r2 
1 

4 



(3 - q; - - CsH^) InaV - 2.80639 + 1.43576q; 



-0.49196 cs^y + 0. 10443 4^ 



+ci ( In aV + 0.86075 - 0.30348 Csw 



+C2 ((g - a) InaV - 33.31690 + 4.29201q; 



+ 1.16439 CsvK + 0.04840 c^j^ 



+ ^am-[-f^p^ + 7^p^) — ^ 
2 Ibn^ 



13 1 \n 2 2 

y + a + - Csvy 1 In a p 
+1.76784 - 2.93576a + 2.56080 c,^,, - 0.92094 



+ci (InaV - 2.06261 - 0.62935c 



sw 



11 



+C2 ( ^ In aV + 0.56785 - a - 4.93253 Csw - 0.19761 



+ ami 



7 2 
- + 2q;+ Csw + C2{- - 2a) 



P^ 16 TT^ 

The contact Green's function is given by 



C"^(p, m) = A'^(p, m)5'(p, m) + S{p, m)A"^{p, m) , 



-A-o, 



where 



(S) 1 

A^'s (p, m) = i-(7MPi. + 7<^Pm) 



+ij(7^p, + 7.P/.)C^ (InaV - 14.27168 - + 0.19740 c.^^) 

Z 16 TT^ Z 

ml, , .Sf^Ci? 1 , ^ . 
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^(6) 1 

+i^(7MP. + luPi.)(^ (InaV - 14.27168 - + 0.19740 c,^) 



(A.31) 



In the MS scheme we have 



.1, ^^'C'i. 

+12(^ + 7^.^/.)^^ 



1 



+i 



P^ 16 TT^ 



,8 p2 31 

-m^^4. 



16 TT^ 



(A.32) 



A.l First unpolarised moment, diagonal, traceless 



In the Ti^^ (i.e. diagonal, traceless) representation of the hypercubic group is 



(3) , . 1 - 1 / ^ ^ \ 



2 " 2 



A 



(A.33) 



In this expression, repeated and v indices are not summed over: 



(3) 1 1 

A^i (p, m; 0, ci, C2) = i- (7/.Pm - 7vPi.) + C2a-{p^p^ - p^py 



+ i^(7MPM-7z.Pz.)f^ 



;- - a) InaV - 17.52702 + 3.79201q; 

3 



-1.72093 csw + 0.35754^,^ + ci (-4.27417+ 1.08793 Csw) 
+C2 (-6.67330 + 4.53710 c,„,. + 0.44621 4^) 



2 16 TT"' 



-3-". 



— m- — „ — — - 4 



p^ 



16 TT^ 



+ a-{p^p^,-p.p.)j^ 



(3 - q; - - Csw) In a^p^ 
-4.58870 + 1.43576 a + 1.14260 Csw - 0.07987 c 



sw 
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+ci (^-^ In aV + 1.88483 + 1.45305 Csw^ 
+C2 (^(^ - a) InaV - 35.68903 + 4.29200 



a 



+0.97769 Csw - 0.04918 c;^ 
+ 10^2^^'^^'^"^"^''^ 16^ 



13 
~6 



a 



i CsSj In a V 



-1.00108 - 2.93576« + 2.62151 c,w - 0.74162 
-ci (inaV - 2.06261 - 0.62935 c^v 
11 



+C2 (^y InaV + 0.38293 -a- 4.95277 c^w - 0.201684^^ 



+ ami- 
2 



16 TT^ 



7 2 

- + 2a+ Csw + C2(- - 2a) 



(A.34) 



The contact Green's function is given by 



C^{p, m) = A^{p, m)S{p, m) + S{p, m)A'^(p, m) , 



.o 



-A-O, 



where 



(A.35) 



(3) 1 



1 a'^C 
+^2^1^^P^i - ^-P-)j^ (1"«V - 14.27168 -^a + 0.19740 c^^) 



1 

2' 



m 



1 



p^ 167r^ 2 

^(3) 1 

A 1 (p, m) = i-{lf^p^^ - luPu) 



-l + a) 



{AM) 



16 TT^ 



(InaV - 14.27168 - -a + 0.19740 c^vk) 



m X 



^^i>7^{liiPlx-lvPv) 



g'Cp 1 



P 



167r2 2 



- -l + a . 



(A.37) 



In the MS scheme we have 



Ams(p, mMs) = i-(7MP/^ - IvPv) 



+'-\{i^P,.-1uPu)y^ 

2 16 TT^ 



- a) In — - 



a 

3 



31 
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-m- „ -— — - 4. 



p. 



167r2 



(A.38) 



A. 8 First polarised moment, off-diagonal (^^v), symmetrised 



For the representation we use the operator 



(6) 1 - 1 / <-> \ 

^]^u = (1 + « ^ Co) - 2 (7^.75^1. + 71^75^^; ^ 

1 - 1/^2 

--mci ^/'CT^v752 [Dy-D^Ji) 



Yl '4^l{<^f^M5{Dx,D,} + a,^^,{Dx,D^})il; . (A.39) 



Repeated hot u indices are not summed over. The amputated Green's function 
is: 



(6) 1 

A""* (p, m; 0, ci, C2) = i^ (7m75Pi/ + 7i/75P/.) 

+ ciai^(7^^75(p^p^ - + caai^ ^ {cTuxlbPuPx + cr^xlbPfiPx) 



.1, x^'C-i. 

+ 12 (7m75Pi^ + 71^75^^) 



- Q;)lnaV 



-19.74374 + 3.79201 a + 0.88956 Cc 



0.49529 ct 



+ci (-5.61603 + 4.10778 c^v^ - 0.26315 4^) 
+C2 (-8.29791 + 4.21724 Csw - 0.49384 4^) 



+ i 



-^75 



167r2 



a 

3 



1 , 

^.t...,o'k I — n..n..] 



' ib TT 



2„2 



- — a + - Csvy) In a p 
2 



3 3 

-4.96628 + 2.43576 a + 0.86287 Csw + 0.06741 c^^ 
+ci ((- - a) In aV - 34.73952 + 3.62534 a 

-0.23227csv^ + 0.017774^) 

/ 5 2 
+C2 -- In aV + 2.39179 + -a + 0.14576 Csw - 0.00452 c; 



3 " 3 



1 



+ "^o m {<y,J.XlbPuPX + CFuXlbPuiPx) 
X+jx,v 



-3.13860 + 2.43576Q; + 0.28852 c^vk - 0.00271 c 



sw 
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+Cl 



-^InaV + 1.53155 + -a 



0.30307cc,^ + 0.01158 ci 



+C2((- - a) InaV - 32.34394 + 3.62534q; 



-0.43312 csv^- - 0.03719 

. 1, ,9^Cf 



16 TT^ 



[-2 + 2q;] 



+ a - ^ Cs,v ) InaV 



+4.77328 - 3.93576a; - 1.36644 c.^k + 1.05432 c 



-'SW 



+ci (^InaV - 2.30755 - - 2.28022 c^vk + 0.07834 c; 



+C2 



- In a V + 0.38901 --a- 2.29812 c^w + 0.21844 4^ 

3 3 



16 



1 + 2q; + - Csp^. 



. 1 m / o o\ 



The contact Green's function is given by 



(ci - C2) 



3(ci + C2) 

8 4 ■ 

a 

3 3 . 



(A.40) 



-irO, 



C (p,m) — A (p,m)S(p,m) + S(p,m)A (p,m), 



(A.41) 



where 



A^i'' (p, m) = i^(7M75Pi. + 7i^75P,.) 



+ i2(7M75Pr. + 7r^75P/.) 



m 1 
'f2 



{lulbPv + lulbPn)]^'- 



„2fi 1 

— ^ (InaV - 14.27168 - -a + 0.19740 c<j,^.) 
g'Cp 1 



167r2 2 

^(fi) 1 

A * (p, m) = i-(7M75Pi. + 7!^75Pm) 
+4(7M75Pt. + lulbPn)^r^ (InaV 

2 16 TT^ 



;-i+«), 



(A.42) 



14.27168 



-a + 0.19740 Csw) 



p2 2^ 



167r2 2 



(A.43) 



In the MS scheme we have 
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.1, ^q'Cf 

I iO TT 



,8 , 31- 

aim— 

-3 ^ //2 9. 



+1 



p 

, 1 m 



-^75 



167r2 



1 



— a 



167r2 



[-2 + 2 a] 



(A.44) 



A. 9 First polarised moment, diagonal, traceless 



For the rj^^ representation we use the operator 



(3). ,1-1/ <-> <-^\ 

(^li = (1 + a m Co) - V 2 (t^xTs^/. - li^lsD^ ) V' 



-^aici 5^^/;^ (^a^xl5{Dx,D^} -a,xl5{Dx,D,})^ . (A.45) 



Repeated /i or u indices are not summed over. The operator vertex is: 



(3) 1 

A""-* (p, m; 0, ci) = i^ (t^xTsPm - T^^TsPi^) 
1 



2 16 vr^ 



- a)lnaV 



-19.92148 + 3.79201 a + 0.99934 c.^^ - 0.60077 4 



sw 



+ci (-15.31376 + 8.54773 Csw - 0.26036 cl^) 



+ ^2 ^^^^^16^ 

+ ai^ ^{(^nXlbPuPx - (^uXlbPvPx) 



1 

3-". 



167r2 



-2.72696 + 2.43576 a + 0.35052 Csw - 0.127617 c^^ 
+ci ((1 - a) In aV - 30.32684 + 4.29201 a 

-0.64293 Cs,^.- 0.005974, 



Im^ 



16 7r2 



[-2 + 2 a] 



1. .g'Cp 



+ ami-(7^75P^ - 7^75^,.) 



167r2 



1 7 

+ a - -CsH-) In aV 
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+3.94429 - 3.93576Q; - 1.41190 Csw + 0.90653 4^ 

+ci (^y In aV - 2.33015 - 2a- 4.68331 Csw - 0.03571 c]^^ 



+ ami 



2 p2 -"^igtt 
The contact Green's function is given by 



-l + 2a + - Csw - -ci 



(A.46) 



C^{p, m) = A'^(p, m)S{p, m) + m)A^(p, m) , 



,0 



■0, 



(A.47) 



where 



+iJ(7/.75PM - 7.75P.)f^ (InaV - 14.27168 -^a + 0.19740 c^^^) 

Z ib TT 



1 

2' 



^(3) 1 

A ^ (p, m) = i-(7M75P/. - li^lsPu) 

+iJ(7M75PM - 7.75P.)f^ (InaV - 14.27168 - + 0.19740 



(A.48) 



m 1 



p z 



167r2 
^'Cf 1 



167r2 2 



;-!+«) 



(A.49) 



In the MS scheme we have 



Ams(p,"^ms) = 12(7^75^/. - 7i^75P^ 



+1o(7m75Pm - lu^5Pu)j^ 
Z iO TT^ 

,.1 (Pm-p') ^ 9'Cf 
2 p^ 16 TT^ 



.8 p2 31 
- — am— — 



, 1 m 



zp ^ 



167r2 



[-2 + 2 a] 



(A.50) 
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